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INTRODUCTION 

Recently the analogue and digital computers are frequently 
used for working out complex automatic control systems, invest i- 

gation of the accident and start up conditions of atomic power plants 
/"A, 5, 67. The use of computers is difficult since equations for the 
subjects under control (in particular^ reactor, heat exchangers) 
are partial differential equations. 

Digital computers permit partial differential equations to be 
solved. Complicated and long programming process of the problem is 
a considerable disadvantage of these machines as well as their 
difficult application to the inve stigation when connecting with 
real control system elements. 

The preparation and treatment of machine equations being simp- 
le, it is not difficult for the analogue computers to connect with 
the real elements of control systems. aerial analogue computers (/W-7 , 
3MU~8 t 3My- 10, MH- 10, MH- 14, MH- 8 , MH-B ) are however designed 
to solve usual differential equations. The method for transformation 
of partial equations with respect to three variables into non-li- 
near differential equations is therefore necessary. 

The simulation method of transient thermal processes in gas- 

cooled reactors on the analogue computers. As is known, transient 
heat transfer in fuel elements of the reactor in a general form can 
be presented by partial equations with respect to three variables: 
the radius of the fuel element, coordinates along the length and time. 
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A met nod for averaging or gus und i’uel eleiu «u Cunhloug~bhe raddn a 
and the reactor length is necessary. to transform this equations in- 
to ordinary differential equations. 

Averaging the teiupex’ature dia trioution of tin fuel element along 
the radius can be realized after introduction of the power dependen- 
ce 

T er (i) = a +&x + ci* + 

( 1 ) 

for the set of equations j 

jjg J otr n ( -r r ) 

dT /* / c c / C ' r / t 

i-Zr _ °£jlIL- ( r -T ) - SfiJL lT r 

9V Mr Sr l'"" ' r) C^TSr W ( 3 ) 

where 

^* 3 ^»®r =mean cross-sectional temperature of fuel element anri 
respectively;, 

7^/7 =fuel container temperature of the gas side; 
fyi =heat flux per unit volume; 

capacity and specific weight of gas and fuel element 
material, respectively; 

n charnel perimeter, cross-sectional area of gas and fuel ele- 
ments ; 

<*V = individual heat transfer coefficient 
G- = gas flow rate. 

For a plane wall, & in width, with _the ^aternal sources anl one- 
sided heat sink with initial condition £»/ TmU^/S' and boundary con- 
ditions ( 0 and -A$jf'=« r (7Z„-Tr) ,it is obtained that 

the thermal resistance of the wall is taken into consideration in tte 
reduced individual heat transfer coefficient 

«», = i/( $ *£c) w 

For rod-type fuel elements of the radius r the reduced indvidual 
heat-transfer coefficient is as follows 

°Cn P = i/(^r r + 

where A is the thermal conductivity of the fuel element materials, 
and the set of equations (2) and 0) becomes 

^ S m ^ O’? ftp /7 f ~T~ — 7” \ 

C„t„v„ ‘ M ''-J 
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/' 


After averaging fuel clement and gas temperatures along the ra- 
dius, the heat-transfer process in the reactor, regenerator and 
cooler is described by partial differential equations with respect 
to two variables (time and the coordinate along the length) 

It is necessary to find the method for determination of para- 
meters along the reactor length. Averaging temperatures along the 
reactor length can be presented by the relation between the moon 
gas temperature and inlet and outlet temperatures 
Trcp - ( Tjp + Ti)/ 2. 

In this case the core must be divided into so many sections along 
the reactor length that within a section the temperature could bo 
considered linear. However, the finite difference method cannot be 
recommended for study of dynamics of an atomic power plant, because 
in this case the sets of equation of the transient heat transfer in 
the heater and reactor occupy a large part of the machine blocks. 

Substitution of the set with distributed parameters by the set 
with lumpdd parameters [^, 5 ] is the simplest and widely used method 
of averaging temperatures along the reactor length. This method is 
admissible at small diversion of the parameters from their nominal 
values. Accuracy of the solution, is uncertain in the case of great- 
ly changing gas flow rate, thermal capacity of the reactor and in- 
let gas temperature. When heat exchangers are presented by a set of 
equations with lumped parameters in the case of great changes of 
variables (start-up, heating, emergency conditions etc.) the tran- 
sient process duration is distorted; the number of peculiarities 
of "the heat transfer process is not taken into consideration. 

The following averaging gas temperature method permits to ob- 
tain an admissible accuracy of simulation with a smaller number 
of operatio na l amplifiers than in the case of dividing the core 
into the sections when the difference method is applied [ 7j • 


Gas-cootdd reactor with rod-type fuel elements. Fig.l shows 
the design diagram of the core with the rod-type fuel elements. 
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The usual assumptions inquired, fox* tstudying transient heat trans- 
fer between gas flow in the channel and the wall of the channel* 
with internal sources of variable heating rat© are following: 

l*Gas subcooling in both bottom and top reflectros is not 
accounted for* 

2*Heat conducted in tho wall along the channel is neglected* 

3* Experimental individual heat-transfo* coefficients of the 
steady-state flow is used for transient processes* 

4*Heat release varies according to a cosine curve along the 
length of the fuel element* 

With thfcse assumptions, transient heat transfer in reactor fuel 
elements is described by the following set of equations: 


CijU, $ =• (cos N - «n ( T„ - Tr) 

c '^ Sr 'If +W 


(2b) 


(3b) 


The channel centre is a zero reading along the length. 

Initial and boundary conditions are: 

l=-*/z Tr(z,~ t/s) = Tap a at Z-o 

2 = e /z Tr (z, e/z) = T, q at t --(? 

Tr(z~ e/z)= ^ 

Equations (2b), (3b) are integrated within the limits from --/z 
to l/z 

tfa ft 

J T„a,i)cU = ~ 




t/i 


’ft a t/i 

- <xn(jT M c^i - JT r c/ 2 ) 


ft 


-ft 

ft 


-ft 


( 6 ) 


( 7 ) 


fyrcfrjt f T r^ - ^ ( fM -/ T r cfz) - Cp <? (T 3 -T„) 

'** v/t)fw>irW =mean gas tempex n a"cure along tJse 
channel length; ^ ( 4 /e) = Tmc P = mean fuel element 

temperature alibng the channel length etc* 

Taking into account for these notations equations (2b) and 
(3b) become of the form: 


d T~mcp _ & ^ , £ +&h. j 


r^cp _ — ) a/- ctF- /r- 7~ ) 

d v Cj„v„ xev** £+&a.) n cfjfrJK*' 
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dTrjp 

ctx 


*F 

‘ptfrcp'Vr 


( 7%cp~ frcp) 


- £& Q. 

K rep Vr 


( Tt Tip) 


(7a) 


where 

v = volume of fuel elements or gas per fuel element; 

P = heat transfer surface of fuel elements; 

Thus for mean temperatures of the gas T rcp and. wall Tncp , 
ordinary non-linear differential equations are obtained. The analo- 
gue computers can be used for their solution. 

Let us apply integration by parts within the limits from ~ p/ 2 
to % to integral of the form J '^T r (Z,Z)dZ to obtain 

the relation between the mean gas temperature T rt p and the inlet 
~f 2p and outlet 7j, gas temperatures. As a result of successive 
integration ,the integral can be represented by the infinite series 

r - <[ k T(r?)d7- + J£l'- 

Trcp-r[ e/ Tr(T,i)dZ-jTr ff ^ 


_zl <l! 

2W dZ 3 S20t dZ* 


l /i 

U/ z 


I 


Infinite series (P' may be presented by the general term of 
the series: 




Z° ) d n T r , .\*+* 


~$T r 


HY 


(as) 


For this kind of series uniform convergence takes place, if ttee 
series converges uniformly to the set z, and the func- 

tion * (at any z and n) forms the monotonous sequence and 
at any z and n is limited 

l&l* * 

B ': the power series converge at any z (the convergence ra- 

d ' .'-I E= <*o) /I - 2/ • 

The function T r (t>Z) is continuous within this series integra- 
tion limits - % .Derivatives along the coordinate z T r Z) 

are limited as the heat flux and mean gas temperature are fiite 
(integral (8)). Therefore infinite series (8) is uniformly conver- 
gent one, and the solution of the problem may be represented by a 
few terms of the series expansion of T, • 

(a) One term of series expansion. 

T ^ri v r 'S't l l v * r>) <*> 
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r 


Hence 


^5 ^ 7*rep ~Tzp 


A set of equations: 

ThC£ 


= r~h ~ ~(i<* f T -r ) 

tt e ' B*zA ' c y s ' ' rc /v 

/y p/y ^ 


(6a) 


^ T. 




IS£ - 


o<F 


' m r<fr<p\ / r 




^3 ~ 2 ^~ rc fi ~ Ftp • 

(bJlVo terms of series expansion 

ft 


(7a) 


2 yt 

T-‘lr.-fy%+ ...l^r 3) . i[(m . (W J 


re s> ~ e 


ae values ((\^) 3 and CQX^m obtained from equation (3b); ^ ^ 


{ (W, 




olt 


Tr',’i(T. r TTj-tftffc-T'-rntfy.&gtysg .]& <g v j 


Hence 


^ l r1 i ’ 8 w^ Trcf T *r) 

'PhA <3*4+: r\-p /Neve.evX ' * . . . 


Ibe set of equations is obtained: 


(6a) 

(7a) 


> * -£U^ - *w 

^ X fop djlf j. fp G & /t - __ Tip A 7j ) ^ /y „T~T ,7“ ) 

jr 3 eftr c r ^v r T re r~ « / ^ ” £ ^^(io) 

(ii) 


°tlz} = (crSK e/f * )N - °< f - fT T ) 
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Hk - (tosx — )Ff- - — (r HZ - To p } 

dl e+2h ' P ' ( 12 ) 

(c) Three terms of series expansion 

T =J-T ^ 

rc/>- ^ '/• ^ g z ^ ^ ( 8d ) 

^e second derivatives of the temperature are determined by equation 

**r r JJlftTM _ l)Tr\ OrfrSr T7? 

TP = CfG- ‘ JT TzJ C/>& ' }2)T 

Sinoe the function T r (i, ijia continuous, has all kinds of partial de- 
rivatives (they are also continuous) and is defined within the 
open range - e /z + ^/z ,the following relation is valid: 
ii Tr_ 

Substitution of integration limits (into equation (8d) and simple 
transformations yield 

T _ / /-r r / / r &Sf s , Cpfo'Vr dp T * 7 

' r*p- 2 * ( 'tp+'d ~J • c P & ty.Q. ' j-p 'rej+ 

i rj P /—' -r‘ ) Cp^rcpr d]rt 1 

+ U L CpS- ( Tf1i Tr */ C/>& ' dz J , H v 


4S L Cpf- 

fhe set of equations: dJhi 

dr 




( t mi~ T ri ) 


ijd. - « 0 *e%^r 7 ‘ 

jf-=- Tr ' f Oi 


£ ' f; J.F 

dl ~ CrfrcpVr 




d Tn ep " _ df-2 ft 

dz C* fa V h ’ jre 1 


/ . # , _ dF / _ , 


Q j"* - 2*4 V# f r ^ £/z I -7 1- /T" t ) 

dr fa* -ft ■ v„ ( CcsS FTdrj d - c„r M i/„ ( r ”> - V 

dl~Ml V* fnath ! 77 - fr — 

dr ‘ fv&ki L e TiTTjd CQZT rm-ty 


Cm fa 
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Analogous sets of equations muy be writcen for 4-,? etc terms of 
series expansion. 

Results of electronic simulation of transint heat-transfir equa- 
tions .The above suggested method permits to describe the transient 
heat— transfer process in apparatuses with various accuracy depending 
on a number of the terms of the series expansion. A larger number 
of terms of the series expansion suggests more accurate mathematic 
description of the process. However, we obtain the set of differenti- 
al equations of a higher order that leads to a sharp incmsase 
of simulating apparatuses ,!£he investigation of transient heat- 
transfer process in a gas-cooled reactor under various external 
perturbations in case of one, two and three terms of the series 
expansion was carried out on the analogue computer of the type 
1 wihtthe aim to choose an optimal number of terms of the 
series expansion. Reactor kinetics have been described by known 
equations (4) taking into account the six groups of the delayed 
neutrons 


dJj. n f r A . r. 

dz ' e* n ■ AlLl 


6 

z 

<</ 

A 




wliere 

n=tfie neutron density 
Ak=the reactivity 

J3 t = the delayed-neutron fraction of the i-group in the total 
number of neutrons 

the decay constant of source of i-group delayed neutron 
Ci = the concentration of nuclei-delayed neutron emitters 
<f*=mean neutron life-time* 

The effective life— time 1*=10 ^sec and zero temperature co- 
efficient for both the fuel and the moderator of the gas-cooled 
reactor were assumed. 

Verification of the method proposed was carried out for 
gas-cooled reactor of an atomic energy plant 50 Mwt in power 
with rod-type fuel elemtns.The coolant was nitrogen*Fuel elements 
composition IK^+BeO, moderator was Be 0, fuel element coatings 

were stainless steel alloys. 
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The estimation of the accuracy of che relation f ^ in 
terms of inlet and outlet temperatures and their derivatives 
along the coordinate z by means of various numbers of tarns of 
the series expansion of T rep at stationary and transient condi- 
tions for the gas-cooled reactor of the atomic energy plant N^= 

50 Mwt was made. 

The estimation of the error of the truncated series can be 
made with sufficient accuracy by the first term of this series 
( 1 ). 

In tire steady state case the accuracy of the series expan- 
sion of the relation T rc/ , in terms of the inlet and outlet 
temperatures and their derivatives, can be estimated in the 
nominal regime. As calculations showed, two terms of the series 
expansion provided sufficient accuracy in the steady-state case 
(the error was 1-5 per cent). 

The transient conditions of a gas-cooled reactor 
50Mwt have heen stodied on the analogue computers of type MH6 -1 
when T rcp was presented in terms of inlet and outlet temperatu- 
res by various number of terms of the series expansion and 
by the finite difference method. 

Analogue cc • .-r MH6 -1 is a block-type arrangement. Ope- 

rational ampli:; iers of this computer have the automatic stabi- 
lization of zero level (the modulation-demodulation channel) 
and, thus, tlie little outlet drift (40-60 mkv for 10 min)The 
frequency responce of the amplifier is linear in inverted 
regime to a frequency of 1000 c/sec. Constant coefficients are 
given by tlie compensating mechod with an error of 0.2 per cent. 
Recording of the results is made on the electronic automatic 
otentiometer;the frame passes the whole scale for 1 sec. The 
.vtic error of the recorder is 0.5 per cent. 

For a gas-cooled reactor on an analogue computer MH6-1 
tue discontinuous and sinusoidal disturbances separate i* each 
parameter have been considered at various numbers of term of 
the series expansion intgansient heat-transfer equation to de- 
termine the limits of validity of this method and the reqdred 
number of series terms. The results obtained are illustraledL in 
Fig.s 2-6. 
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pig, 2 shows the outlet: gas temperature change of the reactor 
with the discontinuous inlet gas temperature drop by 100°. Two and 
three terms of the series expansion give practically the sains values 
flask-ups of the inlet gas temperature. One term of the series ex- 
pansion causes 5-fold excees of the outlet gas temperature flask- 
up. with such disturbances. 

Pigs, 4 and 5 show the change of the neutron power of the re- 
actor and outlet gas temperature at discontinuous decrease of mass 
gas flowrate to G =0.25 0 Q . Flash-ups of the outlet gas temperature 
J practically coincide with two and three terms of series expansi- 
on. The outlet reactor gas temperature flask-up is overextimdred 1.4- 
1.5-fold with one term of the series expansion T rcp . 

Under dropping of the safety rods and constant mass-fhv rate 
of a gas and inlet gas temperature of the reactor, the curves of 
the transient process dor all variables (n,f„ £/ ,,f 3 ,T rc/) ) show a good 
agreement with each other with 1,2,5 terms of the series expansion 
T . The mean metal temperature change in this case is shown in 
Pig. 5. Under the reactivity disturbance, the transient heat-transfer 
process may be described by equation with one term of T rtp expansion 
in terms of inlet and outlet temperatures. 

The following dunamic characteristics of a gas-cooled msactor 
with a complex of perturbations were studied, after researching of 
individual influence of the disturbances aG, aT 2 p , h k to estimate 
the accuracy of representation of the transient heat-transfer 
equations by one, two or three terms of the series expansion: 

1. The deep change of the reactor heat power (ak=-0.005) 

2. Eiscontinuous decrease of the mass-flow rate of a gas 
aG=0.?5G . 

u 0 0 

3. Constant discontinuous temperature drop by 100 C or 50 

increase of inlet gas temperature, 
fne result of simulasiou have shown that expression of T rcp 
in terms of T. and T- by one term of the series expansion under 
such, xomplex disturbances gives essential excess of flask-up of the 
outlet reactor-gas temperature and 1*3-1- 5-fold underestimation 
of the duration of transient process in comparison with the real 
values • 

Two terms of series expansion of T rcp in terms of Tp and 
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and their derivatives along the coordinate z provide sufficient 
accuracy of the solution under reactor transient conditions (er- 
ror is 3-4 per cent). 

The same system with Cue couqjlex perturbation have been 
solved numerically by the finite difference method. The results 
of this solution are the criterion of corect simulation of the 
transient heat-transfer equations with various numbers of tie H Kp 
series, expansion, terms. In the finite difference method, the reac- 
tor fule-element length is broken into a few sect ions. Ehelinear 
temperature profile is assumed within each section. 

The channel length is broken into a few sections, the number 
of which is chosen depending on the required accuracy of the so- 
lution (2»4,6 sections). The results of simulation have shon that 
division of the fuel element length into 4 sections provides 
sufficient accuracy of the solution (error is not more than 1- 
3 per cent). The points of control numerical solution by tte fini- 
te difference method are plotted in Pig. 6 which gives a close 
agreement with the curves of the solutions for 2 and 3 terms of 
the series expansion. 

CONCLUiilOh 

The proposed method of simulation of transient heat-transfer 
processes in a gas-cooled reactor on the analogue computer takes 
into account parameter distribution along the reactor length 
by the expression of the integral-mean gas temperature in terms 
of inlet and outlet temperatures and their derivatives as 
several terms of the series expansion. 

'^he method proposed is verified with deep separate and 
simultaneous disturbances of the gas flow-rate, heat power and 
inlet gas temperature. For a gas-cooled reactor of an atomic 
power plant two terms of series expansion provide sufficient 
accuracy of the dolution (error is not more than 3-4 per cent). 

The method proposed is compared with another simulation me- 
thod - the finite difference method- taken as a reference cne 
and it is shown that the methos presented is more economic than 
the finite difference one (the reactor ia described by 5 diffe- 
rential equations in comparison with ti equations in the fhite 
difference method). 

This method may be used for study of dynamic characteristics 
of a reactor with energy change, heating-up starting-up and reac- 
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tor shut-down, cooling. 

The method may be effectively used for study of the dynamic 
characteristics of a regenerator, cooler and other heat exchangers 
of nuclear power plants. 
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FIGURES 

Fig.l Calculation shceme of the reactor core 

Pig. 2 Outlet jgas temperature variation with jump of inlet tempe- 
rature by 1QO°0. 

Pig. 3. Plot of the neutron power with jump of gas-flow rate by 
G =0.25 G 0 . 

Pig. 4 Plot of the outlet gas temper ature with jump of gas-flow 
rate to G =0.25 0 q 

Pig. 5 Plot of the mean metal temperature with dropping of safety 
rods (Ak=-0.O15). 

Pig. 6 Plot of the outlet gas temperature with simultaneous change 
of the flow-rate and reactivity ( A G =0.75 G Q . Ak=-0.OO5) 
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